Coulomb Correlations and Magnetic Anisotropy in ordered Ll CoPt and FePt alloys 
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We present results of the magneto-crystalline anisotropy energy (MAE) calculations for chemically 
ordered Llo CoPt and FePt alloys taking into account the effects of strong electronic correlations 
and spin-orbit coupling. The local spin density + Hubbard U approximation (LSDA+U) is shown 
to provide a consistent picture of the magnetic ground state properties when intra-atomic Coulomb 
correlations are included for both 3d and 5d elements. Our results demonstrate significant and 
complex contribution of correlation effects to large MAE of these material. 
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Recent progress in fabrication and characterization of 
the magnetic nano-particles and thin films based on the 
ordered LIq CoPt and FePt alloys renewed interest in 
understanding the mechanisms contributing to the large 
magneto-crystalline anisotropy of these materials. While 
the earlier studies were primarily motivated by perma- 
nent magnet applications, current research efforts are fo- 
cused on the use of CoPt and FePt alloys for the high den- 
sity magnetic recording. The large and controllable mag- 
netic anisotropy energy (MAE) is then a crucial property 
to overcome a superparamagnetic limit Q. 

The chemically ordered Ll phase of FePt has large 
uniaxial MAE (« 6.6 * 10 7 erg/cc) which is almost by 
two orders of magnitude higher than their disordered fee 
phase (a dramatic MAE increase with chemical ordering 
is also typical for other alloys of this kind: CoPt and 
FePd). The cubic symmetry is broken in Llo phase by (i) 
stacking of alternate planes of 3d element (Fe or Co) and 
Pt along [001] direction; (ii) by tetragonal distortion due 
to 3d and 5d atomic size mismatch. For the tetragonal 
crystal with uniaxial symmetry, the MAE depends on 
polar 9 magnetization angle |2j as, 



E = K X sin 2 9 + K 2 sin 4 9 



(1) 



where K\ and K 2 are the anisotropy constants. The 
series Eq.Q is rapidly converging and K\ >> K 2 for 
CoPt and FePt LIq alloys H|J The MAE is then com- 
puted as the total energy difference when magnetization 
is oriented along [110] and [001] crystal axes (MAE = 
£[110] - £[001]). 

A significant amount of work has been done to calcu- 
late this energy difference from first principles, employ- 
ing the local spin density approximation (LSDA). For 
the case of CoPt alloy, the fully relativistic Korringa- 
Kohn-Rostoker (KKR) method yields the MAE of 0.058 
meV/f.u. which is very different from results of lin- 
ear muffin-tin orbitals method (LMTO) calculations 
in atomic-sphere- approximation (ASA): 2.29 meV 0, 
1.5 meV U, 2 meV 0, the augmented-spherical-wave 
(ASW) result of 0.97 meV 0, and full-potential LMTO 
(FLMTO) calculation results of 1.05 meV [Hj and 2.2 
meV [llj. As can be seen, there is an apparent and 



significant scatter of the theoretical results, while the ex- 
perimental measurements are known to be quite accurate 
and consistent yielding the MAE of 1 meV [3j . The situa- 
tion is similar for the case of FePt alloy: the LMTO- ASA 
calculations yield 3.4 meV 0], 2.8 meV 0, 3.5 meV ||, 
ASW yieldsja 2.75 meV, and FLMTO 3.9 meV and 
2.7 meV [lfj, and none of these calculations reproduce 
the experimental MAE of 1.3 - 1.4 meV for the bulk 0], 
and 1 — 1.2 meV for the films (when extrapolated 
toT=0|I3). 

The accurate ab-initio calculation of the MAE in itin- 
erant ferromagnet is a very difficult task due to its no- 
torious sensitivity to numerical details [l4|. More im- 
portantly, the LSDA (or generalized gradient approxi- 
mation (GGA)) which is conventionally used in the first- 
principles theory, lacks proper orbital polarization due to 
Coulomb correlation effects 0. In this paper, we wish 
(i) to clarify the ability of the LSDA theory to reproduce 
the experimental MAE for CoPt and FePt alloys, when 
highly accurate full-potential relativistic linearized aug- 
mented plane wave method (FP-LAPW) is used for the 
total energy MAE calculations; (ii) to go beyond LSDA 
and to investigate the role of electron correlations. We ac- 
count for the on-site Coulomb correlation effects by using 
the LSDA+U approach and to show that correlations 
should be included for both "magnetic" (Fe and Co) and 
"non-magnetic" (Pt) sites to describe consistently the 
magnetic ground state properties, such as the MAE and 
spin and orbital magnetic moments. For the first time, 
we quantitatively demonstrate a significant and complex 
character of the intra-atomic Coulomb repulsion contri- 
bution to the MAE of itinerant ferromagnet with strongly 
magnetic 3d and non-magnetic 5d-elements. 

We start with the conventional LSDA band theoretical 
method together with the relativistic FP-LAPW method 
and apply them to perform total energy electronic 
and magnetic structure calculations for the magnetiza- 
tion fixed along each of [110] and [001] axes, respectively 
and the MAE. The special fc-points method is used for 
Brillouin Zone integration with the Gaussian smearing 
of 1 mRy for fc-points weighting 0] • For convergence of 
the total energy differences within desired accuracy (bet- 
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TABLE I: Magnetic Anisotropy Energy (meV/f.u.) calcu- 
lated using FP-LAPW method within the LSDA theory. 
Ref. [91 Ref. [Ill Ref. flOl Present 



CoPt 0.97 
Exp. 0.82 (T : 



2.2 1.05 
293 K) [3J; 1.00 (T - 



1.03 

K) [13J 



FePt 

Exp. 



2.75 3.9 
0.7 - 1.2 (T = 



2.73 
293 K) [4, 12} 



1 — 1.4 (T = K) [13] 



ter than a few neV), about 11000 fc-points are used (the 
MAE as a function of Appoints number is calculated to be 
very similar to the results of Ref. 10]). The experimental 
values for the lattice parameters (a — 7.19 a.u., c = 7.01 
a.u. for CoPt, and a = 7.30 a.u., c = 7.15 a.u. for FePt) 
are used 0|- The calculated MAE is shown in Table I. 
in comparison with recent total energy calculations and 
experimental results. The present FP-LAPW results are 
in very good quantitative agreement with ASW results of 
Oppeneer [t| and FLMTO results of Ravindran and 
disagree substantially with FLMTO results of Galanakis 
[l^ for the reason which is unclear to us. For CoPt alloy 
the MAE is calculated in a very good agreement with 
experimental data, while for FePt our LSDA results (to- 
gether with those of Refs. overestimate value of 
MAE by a factor of two. 

Very recently, the electron-electron interaction was 
shown to play an important role for the MAE in itiner- 
ant d ^3 and /-electron magnetic materials. Here, 
we use the LSDA+U method combined with relativis- 
ts FP-LAPW basis [HII3 

to account for the intra-site 
Coulomb repulsion U. Minimization of the LSDA+U to- 
tal energy functional with SOC treated self-consistently 
[l7l ] generates not only the ground state energy but also 
one-electron energies and states providing the orbital 
contribution to the magnetic moment. The basic dif- 
ference of LSDA+U calculations from the LSDA is its 
explicit dependence on on-site spin and orbitally resolved 
occupation matrices. The LSDA+U creates in addition 
to spin-only dependent LSDA potential, the spin and 
orbitally dependent on-site "+Z7" potential which pro- 
duces the orbital polarization [2l|. Since the LSDA+U 
method is rarely applied to metals, the appropriate val- 
ues of intra-atomic repulsion U for Fe, Co and Pt atoms 
in metals are not known precisely as they are strongly 
affected by the metallic screening. Here we choose U^d 
from the range of 1-2 eV, and U§d from the range of 
0-1 eV (according to the experimental values for pure 
metals |22J) in order to reproduce experimental values 
of the MAE. For the exchange J we use the values of 
J Co = 0.911 eV, J Fe = 0.844 eV, J Pt = 0.544 eV which 
are obtained as a result of constrained LSDA calculations 
[23| . and are close to their atomic values. 

First, we discuss the FePt alloy. The MAE as a func- 
tion of U Fe is shown in Fig. 1(a). When U is included 
on Fe-site only and varied in the interval of 1-2 eV, the 
MAE is decreasing from its LSDA value. The experi- 



2.5 1 



> 


2 







E 






1.5 


llT 




< 






1 




0.5 



• •+U(Pt) = 

-i 1 — ■ B+U(Pt) = 0.54 eV 

LSDA ♦+U(Pt) = 0.82 eV 



Exp. MAEI 



a) 



1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 

Intra-atomic Coulomb U, eV 



• •+U(Pt) = 

+U(Pt) = 0.27 eV 
+U(Pt) = 0.54 eV 




0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 

Intra-atomic Coulomb U, eV 



FIG. 1: (color) The MAE vs intra-atomic Coulomb repulsion 
(U) on 3d-site for (a) FePt alloy; (b) CoPt alloy for different 
values of U on Pt-site. Note, that we use the bulk experimen- 
tal MAE values 0,3 extrapolated to T = K. 



mental value of the MAE can not be reached without 
unreasonable increase of U. Further meaningful reduc- 
tion of the MAE can be achieved by including U on the 
Pt-site (see, Fig. 1(a)) and for the value of Up t around 
0.544 eV (= J Pt ) the MAE observed for the bulk FePt 
can be reproduced. We note that the correct value of the 
MAE can be obtained in the region of U Fe ~ 1.5-1.6 eV, 
U Pt « 0.5-0.6 eV. 

We choose the {U Fe = 1.516 eV, U Pt = 0.544 eV} as 
a representative material specific parameters set to an- 
alyze the spin and orbital ground state properties. The 
calculated spin M s , orbital Mi and total magnetic mo- 
ments for the magnetization directed along [001] axis are 
shown in Table II. The LSDA+U agrees with observed 
total magnetization/f.u. M tot , yields a slight increase 
of Mf e from its LSDA value of 2.88 /ip and substan- 
tially increases the Mf e from 0.0657 [ip. The change 
of M s and Mi of Pt-site is very small as compared with 
LSDA. The effect of U on the Mi is seen to be similar 
to the "orbital polarization correction" of Brooks et al. 
(LSDA+OP) Hg (see, Table II). This ad-hoc LSDA+OP 
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TABLE II: Spin (M s ), Orbital (Mi) magnetic moments for 
3d and 5d atoms, and Total Magnetic Moment (M tot ) per 
formula unit (hb)', Magnetic Anisotropy Energy (meV/f.u.) 
as results of LSDA+U calculations for FePt and CoPt alloys. 
The LSDA+U calculated M a and Mi values for bcc-Fe and 
hep-Co obtained with the same values for Uzd and Jsd as for 



and Co-atoms in FePt and CoPt alloys. 


FePt[001] 


MAE M tot 


Atom 


M 

1V1 S 


Nil 


LSDA+U 


1.3 3.47 


Fe 


O.UU 


n 1 1 a 


{U Fe = 1-52 eV, U Pt = 0.54 eV} 


Pt 


0.34 


0.048 


LSDA+OP 


2.9 3.36 


Fe 


9 8Q 


n 1 1 n 






Pt 


U.oO 


n ri/is 

U.U4o 


Exp. 


1.3 3.4 








CoPt[001] 


MAE M tot 


Atom 


M 


M, 


LSDA+U 


1.0 2.55 


Co 


1.93 


0.253 


{U Co = 1.70 eV, U Pt = 0.54 eV} 


Pt 


0.38 


0.065 


LSDA+OP 


1.64 2.37 


Co 


1.80 


0.161 






Pt 


0.39 


0.062 


Exp. 


1.0 2.4 


Co 


1.9-2.1 


0.284 






Pt 


0.36 


0.090 


bcc-Fe[001] 






M a 


Mi 


LSDA+U 


U Fs = 1.52 eV 


2.234 


0.085 


LSDA+OP [14] 






2.193 


0.078 


Exp. 






2.13 


0.08 


hcp-Co[0001] 






M s 


Mi 


LSDA+U 


Uco = 1.70 eV 


1.631 


0.153 


LSDA+OP [14] 






1.591 


0.123 


Exp. 






1.52 


0.14 



correction is known to improve the Mi in a number of el- 
and /-electron compounds, but it does not reproduce the 
MAE for FePt alloy. 

We now describe CoPt alloy. The MAE as a function 
of Uco is shown in Fig. 1(b). When U is included on Co- 
site only and varied in the interval of 1-2 eV, the MAE 
increases from its LSDA value and deviates from the ex- 
periment. As in the case of FePt alloy, we need to include 
U on the Pt-site. The increase of U pt yields the decrease 
of the MAE and its correct value is found in the region 
of {Uco « 1.4-1.7 eV, U Pt ~ 0.3-0.6 eV}. We choose the 
set {Uco — 1-698 eV, Up t — 0.544 eV} as a representa- 
tive (note, that we use the same value of U pt as for the 
FePt alloy, indicating its transferability). The calculated 
ground state M s , Mi and M tot for the magnetization di- 
rected along [001] axis are shown in Table II. The M tot is 
calculated in a good quantitative agreement with results 
of LSDA+OP and experimental data ||. There is only 
small increase in M s c ° from LSDA value of 1.81 [1b while 
the Mf° is substantially enhanced from its LSDA value 
of 0.093 [i B (Table II). The change of M s and Mi for 
Pt-site is very small. The Mf° increase in LSDA+U cal- 
culations is substantially more pronounced than obtained 
in LSDA+OP approximation [lOj. 

Recently, the element-specific M s and Mi for CoPt 
were measured by Grange et dl. u sing the X-ray Mag- 
netic Circular Dichroism (XMCD) 0. The LSDA+U 
calculated values for M s and Mi for the Co-site are in 



a good quantitative agreement with XMCD results [25j. 
For the Pt-site, the agreement is not as good. Probably 
it is caused by the use of the atomic-like sum rules to 
extract M s and Mi from XMCD spectra. This proce- 
dure is not reliable for Pt due to substantially itinerant 
character of Pt 5d-electrons. The LSDA+U calculations 
reproduce consistently the MAE, total, spin and orbital 
ground state magnetic moments for CoPt alloy (Table 
II). Both LSDA and LSDA+OP are only partially suc- 
cessful: LSDA yields correct value of MAE but fails for 
Mf° and LSDA+OP improves somewhat Mi but does 
not reproduce the MAE. 

To evaluate further the consistency of the LDSA+U re- 
sults, we performed the LSDA+U calculations for the ele- 
mental 3<i-ferromagnet bcc-Fe and hep-Co with the same 
values of the U^d and J^d as for Fe- and Co-atoms in FePt 
and CoPt alloys (see, Table II). It is seen that without 
any further adjustments of parameters the LSDA+U pro- 
vides very reasonable results for the orbital magnetiza- 
tion in elemental 3<i ferromagnet. These results demon- 
strate that the on-site Coulomb interaction parameters 
Usd are well transferable in the transitional d-metal sys- 
tems. 

As for the choice of the Coulomb-L 7 for Pt, the chal- 
lenge lies in correcting the LSDA orbital polarization 
without harming the exchange splitting which is ex- 
pected to be well accounted in the LSDA. The choice 
of Upt =0.5-0.6 eV looks then quite reasonable. Indeed, 
for the Pt-5ei states having almost equal on-site occupa- 
tions, the choice oiUp t ~ Jpt corresponds to an effective 
Stoner exchange Ilsda+u ~ Ilsda H3 preserving the 
LSDA spin polarization. This allows to ensure that the 
LSDA+U correction contributes entirely to the orbital 
polarization. The LSDA+U method, while proposed to 
deal with the problems specific for the localized states, 
in fact is not limited by this case. This method can be 
used as soon as on-site Coulomb correlation in the form 
of the Hubbard model is physically meaningful. The 
above comparison with available experimental data for 
the MAE and spin/orbital magnetic moments, and phys- 
ically reasonable choice of parameters justify the use of 
on-site LSDA+U correction T6j for the Pt-5<i states as 
the way to correct on-site orbital polarization. 

We now discuss the relation between the MAE and 
the anisotropy of the orbital magnetic moment AMi(= 
Mi || [110] - Mi || [001]). Bruno H| showed that in 
the limit of strong exchange splitting A ex >> SOC, the 
MAE is proportional to AM;. This model predicts the 
positive MAE of 0.2 meV for CoPt and negative MAE of 
-2.1 meV for FePt, in disagreement with the total energy 
calculations and experiment (Table II). 

A more general form for MAE was given in Ref . : 

MAE w -|A(M/ - Mj) + AEt [T |- "spin- flip" ] (2) 
where £ is the SOC constant. The 1st term (AEl) is the 
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Tim-spins contribution due to the orbital moment L, 
and the 2nd term couples the ||-spins and is related to 
the spin magnetic dipole moment T. In the limit of SOC 
<< A ex , the AEt ~ — 3£ 2 / A ex {AQ zz ] is proportional 
to the difference of quadrupole moments Q zz for z — 
[110], [001]. Note, that this AE T form is valid for Fc 
and Co since their SOC (0.07-0.08 eV) << A ex (3-4 eV) 
and can not be used for Pt which has the SOC (0.6 eV) 
> A ex (0.2 eV) [H|. 

The Eq.@ gives for the CoPt alloy A£?£°=1.4 meV, 
AE$°=-0.7 meV and AJ7£*=1.9 meV, and for FePt al- 
loy, A££ e =-0.1 meV, A£f e =1.0 meV, and AEp=1.2 
meV. Here, large AE^ e ' Co contributions to the MAE 
naturally originate from the difference in the inter- 
plane {xz,yz} 3d — 5d and in-plane {xy} 3d — 3d hy- 
bridization. Without AE P * contribution, the total MAE 
of 2.7 meV (CoPt) and 2.1 meV (FePt) can be es- 
timated, exceeding substantially the experimental val- 
ues (cf., Table II). We can only roughly estimate that 
AEp ~ -AQ ZZ (= 0.12 (CoPt),0.09 (FePt) provides ad- 
ditional negative MAE contributions which are expected 
to reduce a total MAE towards the experimental data. 
Thus, due to the strong Pt-SOC, neither of commonly 
used MAE parameterizations [ifl based on SOC- 
perturbation theory expansions is valid on the quantita- 
tive and do not provide a substitute for the total energy 
MAE calculations. 

Still, it is of interest to apply the Eq. (J2J to analyze 
qualitatively the origin of MAE dependence on U shown 
in Fig.l. This analysis shows that LSDA+U MAE vs U 
dependence is qualitatively consistent with the Eq.J2J, 
and the change in AE% * ~ A(Af/ - M}) pt contributes 
substantially to the MAE variations with the U. In par- 
ticular, we find the "coupling" between C/ 3( j and AE pt 
which originates from strong 3d — 5c? hybridization, so 
that tiny L^d-induced changes in Pt-A(M/ - Mj) pro- 
duce substantial MAE change due to the strong Pt-SOC. 
It also explains surprisingly strong MAE dependence on 
Upt (cf., Fig.l), as a variation of Upt causes tiny Pt- 
A(M} — Ml) change (~ 10~ 3 lib) which in turn changes 
the MAE substantially due to the strong Pt-SOC. 

To summarize, accounting for on-site Coulomb cor- 
relations beyond what is included in LSDA, and us- 
ing the LSDA+U method in a very general implemen- 
tation including SOC we have provided a microscopic 
theory of the ground state magnetic properties in LIq 
FePt and CoPt alloys. It is shown by comparison with 
the experiment that LSDA+U method is capable of de- 
scribing quantitatively the MAE and orbital magneti- 
zation in these alloys with physically reasonable choice 
of Coulomb-?/ parameters. Using the SOC-perturbation 
theory model we provide interpretation of our numeri- 
cal results. These results are believed to be important 
for quantitative microscopic understanding of the large 
MAE in these material, and will assist in the develop- 



ment of the next generation magnetic recording devices 
operating above Tbit/in 2 recording densities. 
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